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Abstract — Optimal quantum stabilizer codes of distance 3 are 
explicitly constructed for all lengths except for the following four 
families of lengths 8/m — {1,2} and f,n+2 — {2,3} with fm = 
^ 3^"*" and m > 2 being integer, for which our codes are of the 
best parameters known and are only one logical qubit less than 
the quantum Hamming bound. The optimality of our codes is 
ensured by saturating either the quantum Hamming bound or 
a stronger bound for three families of lengths 8/m + {1,2} and 
/m+2 ^ 1 with m > 1 derived from the linear programming 
bound. For the lengths less than 128 three previously unknown 
codes [[36,29,3]], [[37,30,3]] and [[81,73,3]] have been found. 

Index Terms — quantum error correction, 1-error correcting 
stabilizer codes, quantum Hamming bound, linear programming 
bound, optimal codes 

I. Introduction 

Quantum error-correcting codes [2], [9], [13], [15] provide 
us an active way of protecting our precious quantum data from 
quantum noises and play essential roles in various quantum 
informational processes. Simply speaking, a QECC is just 
a subspace that corrects certain type of errors. When the 
subspace is specified by the joint +1 eigenspace of a group of 
commuting multilocal Pauli operators, i.e., direct products of 
local Pauli operators, the codes are called as stabilizer codes 
[3], [4], [5]. We consider only binary codes here. As usual 
we shall denote by [[n, fc, d]] a stabilizer code of length n 
and distance d, i.e., correcting up to J -qubit errors, that 
encodes k logical qubits. 

One fundamental task is to construct optimal codes, e.g., 
codes with largest possible k with fixed n and d. In the case 
of d = 2 all optimal stabilizer codes are known. In the simplest 
nontrivial case d = 3, despite many efforts to construct optimal 
stabilizer codes, a systematic construction for all lengths has 
not been achieved yet. Known results include Gottesman's 
optimal codes family [6] of lengths 2™ with m > 3 which 
has been generalized for even lengths [10] by using Steane's 
enlargement construction [16] with some codes being optimal 
and some are suboptimal, i.e., one logical qubit less than the 
quantum Hamming bound. 

A code of distance d is degenerate if there are harmless 
undetectable errors acting on less than d qubits, i.e., errors can 
not be detected but do not affect the encoded quantum data. 
If all errors acting on less than d qubits can be detected the 
codes are non-degenerate or pure. For a pure code of distance 
3 all errors happened on up to 2 qubits can be detected. The 
quantum Hamming bound, e.g.. 



k>SH^ [log2(3n+l)] 
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TABLE I 

A SUMMARY OF THE EXISTENCE OF OPTIMAL STABILIZER CODES 
[[n, fc, 3]] OF ALL LENGTHS. ALL THESE CODES SATURATE THE QUANTUM 
HAMMING BOUND EXCEPT FOR 7 FAMILIES OF LENGTHS WITH 3 OF THEM 

(LABELED BY I) BEING NEVERTHELESS OPTIMAL AND 4 OF THEM 
(LABELED BY u) HAVING THE BEST PARAMETERS KNOWN, 
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7,9, 10 


6 
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11 < n < 17,21 


6 




6 




18,19, '20 


7 




6 




22 < n < 35, "36, "37,40 


7 




7 




"38, "39, '41, '42 


8 




7 




43 < n < "81,85 


8 




8 




"82, "83, '84 


9 




8 




86 < n < 128 


9 




9 




+ 1 < n < 8/m - 3 (m > 2) 


2m ^ 


-3 


2m 4- 


3 


P(8/m) (m > 1) 


2m ^ 


-3 


2m 4- 
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8/m - ("1, "2} (m > 2) 


2m ^ 


-4 


2m 4- 
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8/m + {'l, '2} (m> 1) 


2m ^ 


-4 


2m 4- 
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8/m + 3 < n < /„+2 - 4 (m > 2) 


2m ^ 


-4 


2m 4- 
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f /m+2 {m > 0) 


2m ^ 


-4 


2m 4- 
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/m+2 - ("2, "3} (m > 2) 


2m + 5 


2m 4- 


4 


'(/m+2-1) (m>l) 


2m ^ 


-5 


2m 4- 


4 



for a stabilizer code [[n, fc, 3]], is introduced initially for 
the non-degenerate codes, has been proved to be valid for 
degenerate codes of distance 3 and 5 [5] and of a large enough 
length [1] via linear programming (LP) bound [4], [11], 

The most comprehensive list of stabilizer codes of distance 
3 up to 128 qubits is presented in the public code table 
maintained by Grassl [8] and the parameters of the optimal 
codes are summarized in the upper half of Table I, where 
the optimal codes of length labeled by I do not saturate the 
quantum Hamming bound and the optimal code of length 
n = Q (labeled by jJ) is degenerate. Those optimal codes of 
three lengths labeled with a are previously unknown and will 
be constructed here. 

In this paper we shall construct explicitly the optimal 
stabilizer codes of all lengths and a summary is given in 
the lower half of Table I, in which the optimal codes of 
lengths labeled by p are already known. For simplicity we 
have denoted /,„ = ^ (m > 1), i.e.. 



(1,5,21,85,341,...,/, = 4/,_i + l,...}. 



(2) 



All our codes either i) saturate the quantum Hamming bound 
(of unlabeled lengths) or ii) saturate the LP bound that has 
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been worked out analytically for three special families of 
lengths (labeled by I) or iii) have the best known parameters 
and are only 1 logical qubit less than the quantum Hamming 
bound (of lengths labeled by u). 

After the introduction of some notations and known results 
essential to our construction in Sec. II, we shall present the 
general construction for the optimal code of a length n > 38 
in Sec. Ill and then work out the LP bound for 3 families of 
lengths in Sec. IV to ensure the optimality of some families 
of our codes. In Sec. V we construct explicitly all the pure 
optimal codes of lengths n < 37 case by case, which are 
essential for our general construction. 

II. Notations and known results 

Our construction is based on two families pure codes 
and Gottesman's stabilizer pasting [7] to build new codes 
from old pure codes. As usual we denote by X, Y, Z three 
Pauli operators and by / the identity operator Furthermore 
we denote X{n) = XiX2..-Xn with Xi being the Pauli 
operator X acting nontrivially on the i-th qubit only and 
similar expressions for Y{n)^ Z{n), and I{n). For simplicity 
we shall denote by [n, s] the stabilizer of a pure stabilizer code 
[[?!, n — s, 3]] while simply by [n] the stabilizer an optimal pure 
code of length n, e.g., [5] stands for the perfect code [[5, 1, 3]] 
whose stabilizer reads 
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where a juxtaposition of some Pauli operators in the same row 
means their direct product. 

Codes family [2™] (to > 3). The first codes family is 
Gottesman's family of optimal codes [[2™, 2™ — m — 2,3]] 
with TO > 3 that saturate the quantum Hamming bound 
[6]. By construction, these codes are non-degenerate and two 
observables X(2™) = Xi . . . and Z(2") = Zi . . . 
are generators of the stabilizer For simplicity we denote by 
[2™] a set of m + 2 generators of the stabilizer of Gottesman's 
code with the first two generators being X(2™) and Z(2"'). 

An explicit construction of the remaining m generators 
are given by the check matrix [H„i\A„iH,n] where H,n = 
[co,ci, . . . ,€2^-1] with the (fc + l)-th column Ck being the 
binary vector representing integer k {k = 0, 1, . . . , 2"' — 1) 
and Am is any invertible and fixed point free m x m matrix, 
i.e., AmS ^ and A„iS ^ s for all s G F™. As an example 
the unique code [2'^] has a stabilizer generated by 

XXXXXXXX 

ZZZZZZZZ 

I Z I Z Y X Y X . (4) 

I Z X Y I Z X Y 

I Y Z X Z X I Y 

Codes family [8 • m] (to > 3). The second family of codes 
are of parameters [[8to, 8to — Im — 5, 3]] with l,n = \^og2 to] 
that are constructed in Ref.[10]. One crucial property of this 
family is that they are stabilized by the all X and all Z 



TABLE II 

Some examples from codes family [8 • m]. 



I{2^) X{2^) Y{2^) I{2^) X{2^) Y{2^) Z{2^) 
1(2:*) Y{2^) Z{2-^) 1(2*) y(2») ^(2^) X(2-^) 
[8-3] = [[24,17, 3]] [8 ■ 4] = [[32,25,3]] 

\2^] [23] [23] \2^] [23] ^] 
7(23) Z(23) y(23) X(23) y(23) X(23) 
X(23) y(23) 7(23) Z(23) X(23) y(23) 
Z(23) X(23) Z(23) X(23) 7(23) y(2^) 

[8 ■ 6] = [[48, 40, 3]] 



observables X{8m) and Z{8m). Here we shall provide a 
different construction based on Gottesman's codes family. 

We divide 8to qubits into m blocks of 8 -qubit. First 5 
stabilizer of the code are [2^]®™ whose first two generators 
are X{8m) and Z(8to). In the case of m = 3, 4 the codes are 
defined in Table II. In the case of 7ti > 5 so that l„i > 3, 
the remaining /,„ generators of the stabilizer are obtained 
from Gottesman's code [2'™] by at first removing the first 
two generators and then removing arbitrary 2'™ — m qubits 
and finally replacing each single-qubit Pauli operators X, Y, 
and Z in the remaining stabilizers with corresponding 8-qubit 
operators X{2^),Y{2^), and Z{2^) respectively. In Table II 
we also present an example in the case of to ~ 6. 

Obviously all + 5 generators defined above are com- 
muting with each other. Because of the first 5 generators of 
the stabilizer any 2-errors in the same 8-qubit block can be 
detected. For any 2 errors in two different 8-qubit blocks, the 
last l,n generators together with the first 2 generators defines 
a subcode of Gottesman's code [2''"] and therefore detects all 
2 errors in different blocks. Thus all 2-error can be detected 
so that we have constructed a pure 1 -error-correcting code of 
length 8to. 

We shall abuse the notation slightly to denote all the codes 
of this family by [8 ■ m\ though some of them are not optimal. 
In fact when fr+i + 1 < m < 2'^''+^ and ifl^ < m < 2^'' 
with r > 1 the code [8 • to] is optimal since Im + 5 — sh in 
these cases. Otherwise the code is suboptimal, i.e., Im + ^ = 

SH + l. 

Stabilizer pasting (Gottesman [7]): Given two non- 
degenerate stabilizer codes [712,52] = (<S'i, <5'2, • ■ • , <5's2) ^nd 
[ni, si] = (Ti, T2, . . . , Tsi) of distance 3, if two observables 
X{n2) and ^(712) belong to [7i2,S2], say, 5*1 = X{n2) 
and S2 = Z{n2), then the stabilizer defined in Table III 
defines a non-degenerate stabilizer code [77,2 + 77,1,5] with 
s = max{s2, si + 2}, denoted as [712, S2] [> [tii, si]. 

As the first example of stabilizer pasting we can obtain an 
optimal code [13] = [[13, 7, 3]] by pasting the optimal code 
[2"^] of length 712 = 8 and §2 = 5 stabilizers with the perfect 
code [5], i.e,, 77i = 5 and si = 4. The resulting code is of 
length 77i + 7i2 = 13 with si + 2 = 6>S2 = 5 stabilizers. 

If there is a third pure code [773, S3] with X^n^i) and Z{n^) 
belonging to its stabilizer then the stabilizer pasting results in 
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TABLE III 

The stabilizer for the code obtained from pasting. 



X(n2) 


7(ni) 


X{n2) 


/(m) 


Z{n2) 


7{ni) 


Z{n2) 


7(ni) 


S3 


Ti 


S3 


Ti 


Si 


72 


Si 


T2 






or : 




Ss2 




Ssi+2 






T32-I 


Sbi+3 


/(m) 


I{n2) 




Ss2 


/(ni) 



a pure code 

[ni + 712 + n-3, s] = ["-3, S3] > ["2, S2] > [ni, si] (5) 

with s = niax{s3,S2 + 2,si + 4}, which can be further 
pasted with another code and so on. As the second example 
the perfect code [[/,„,./«! — 2m, 3]] with /„i = "Ht^ ^'^'^ 
TO > 3 can be constructed by pasting Gottesman's codes [2^'] 
(Z = 2, 3, . . . , m) with the pure perfect 5-qubit code [7], [4], 

iU = [22('"-i)] > [22("-2)] >...[> [2^] > [5]. (6) 

As the last example the optimal stabilizer code of length 8fm 
(to > 2) can be constructed by pasting Gottesman's codes 
[22'+i] (/ = 1,3,...,to) [4] 

[8/™] = [22™+i] o [22—1] ^ _ ^ ^ ^ p3^^ 
These lengths are labeled by p in Table I. 

III. General construction 

Our main tool is the pasting of codes to produce new codes 
from old ones and only pure codes can be used in the pasting. 
Since the optimal stabilizer code for n = 6 is degenerate the 
optimality does not ensure the pureness. Although from the 
upper half of Table I we know the optimal codes exist for 
n < 37, we have to check case by case that pure optimal 
codes also exist, which are essential to our construction. 

Lemma 1 Non-degenerate optimal 1 -error correcting codes 
of lengths 10 < ri < 17 and 30 < ri < 37 exist. 

Proof: By a direct application of the stabilizer pasting 
to two optimal codes we obtain a previously unknown pure 
optimal code [[37,30,3]] whose stabilizer reads 

[37] = [2^] [> [5]. (8) 

Also it is not difficult to check that the optimal stabilizer 
code [[17, 11,3]] found in Ref.[4] by a random search is non- 



TABLE IV 

The stabilizers of the pure optimal codes [[n, n — s, 3]] for 

n < 37 AND n 7^ 6. ALL THE 2-ERROR-DETECTING BLOCKS SUCH AS 
[28, 7]2 ARE CONSTRUCTED IN SEC. V EXPLICITLY. 
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s 


Stabilizer 


n 


s 


Stabilizer 


10 


6 


Table VII 
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[4,4]i l> [111 

L ' J -L ^ L J 


11 
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[10, 6] 1 t> [l]i 


7 
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[6,6]i > [l]i 


12 


6 


[10, 6]2 > [2,4]2 


8 


5 


[2^] 


13 


6 


[10, 6]2 > [3,4]2 


9 


6 


[6,6)2 l> [3, 4)2 


14 


6 


[10, 6] 1 > [4,4]i 


18 


7 


[10] > [23] 


15 


6 


[10] > [5] 


19 


7 


[18, 7] 1 > [l]i 


16 


6 


[2^] 


20 


7 


[18, 7] 2 > [2, 4] 2 


17 


6 


Eq.(9) (Ref.[4]) 


21 


6 


[21 > [5] 


30 


7 


[28,7]2 > [2,4]2 


22 


7 


[18, 7] 1 > [4,4]i 


31 


7 


[28,7]2 t> [3,4]2 


23 


7 


[18, 7] 2 > [5, 5] 2 


32 


7 


[2^] 


24 


7 


[8-3] 


33 


7 


[28,7]2 > [5,5]2 


25 


7 


[18, 7] 1 > [7,5]i 


34 


7 


[26,7]2 > [7,5]i l> [l]i 


26 


7 


[18, 7] 2 > [7,5]i > [l]i 


35 


7 


[28, 7] 1 > [7,5]i 


27 


7 


[18, 7] 1 > [2^] > [l]i 


"36 


7 


[28,7]2 t> [7,5]i l> [l]i 


28 


7 


[20, 7] 2 > [7,5]i > [l]i 


"37 


7 


[32] > [5] 


29 


7 


[8 • 3] > [5] 



degenerate, whose stabilizer reads 

/ I X I ZYZYXXZYI I X XY 
I IZXIZIYYYXXZYYXX 
I X I IXZXZYYYIYXZIY 
IZIZZIYXYXZYZXZZX. (9) 
XI I Z Y I IXZZYXYZIYX 
Z I IXYYYIYIYXIXXZY 

[17] = [[17,11,3]] 

Obviously pure optimal codes of lengths 16 and 32 exist. We 
shall postpone the explicit constructions of the pure optimal 
codes of remaining lengths to Sec. V where the pasting of 
stabilizers is generalized to the pasting of noncommuting sets 
of generators. In fact all the pure optimal codes of lengths 
5 < n < 37 with n ^ 6 are summarized in Table IV. It 
is worthy of noting that there is another previously unknown 
optimal stabilizer code [36] = [[36,29,3]], whose stabilizer is 
explicitly given in Table VI. ■ 

Lemma 1 ensures that there exist [17 — P] and [37 — P] for 
< /3 < 7, i.e, optimal pure codes of those lengths exist and 
have 6 and 7 generators respectively. For 71 > 38 we have the 
following general construction: 

Theorem 2 For a given length n > 38 if a) 8fm — 2 < n < 
fm+2 ~ 4 f recalling that = ^ ) for some m > 2 then 
we denote fm+2 — 4 — n = 8a + P with a > and < /3 < 7. 
The stabilizer 

[8 • (22™-i - a)] O [2^™] > [22™-2] o . . . o [2^] [> [17-/?] (10) 

defines a non-degenerate code [[n, n— 2to— 4, 3]]. When m = 2 
the stabilizer is generated by [8 • (8 — a)] [> [17 — P]. b) If 
fm+2 ~ 3 < n < 8fm+i — 'i for some to > 2 then we denote 
8/m+i — 3 — n = 8a + /3 with a > and < /3 < 7. The 
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Stabilizer 

[8- (22™-a)] [> [2^™+!] o [22"- 1] o . . .o [2^] [> [37-/3] (11) 

defines a non-degenerate code [[n, n— 2m— 5, 3]]. When m ~ 2 
the stabilizer is generated by [8 • (16 — a)] > [37 — /?]. 

Proof: At first from Lemma 1 and the constructions 
of two codes families [8 • k] and [2^^] it is clear that all 
the stabilizer codes involved in Eq.lfTOli or Eq.lfTTTi are non- 
degenerate. Secondly by construction two families of codes 
[8-fc] and [2*^] are stabilized by all X and all Z Pauli operators. 
As a result the stabilizer pasting can be applied from right to 
left so that Eq.lfTOli and Eq.(fTTTi define pure stabilizer codes of 
distance 3. 

Now we evaluate the parameters of the codes. It is easy 
to see from the definition of a and /? and identity fm+2 = 
22m+2 + 2^™ + . . . + 2-^ + 5 that the length of the resulting 
codes are exactly n. Recalling that the codes [8 • k] and [2*^] 
have Ik = [log fc] + 5 and fc + 2 stabilizers respectively while 
the codes [17 — /3] and [37-/3] have at most 6 and 7 stabilizers 
respectively. Since a > we have [log(22™-° — q)] < 2m — a 
for a = 0, 1, the stabilizers in Eq.(fToli and Eq. lfTTT i have 2m+4 
and 2m + 5 generators respectively. ■ 

Let us look at some examples. If n = 38 then m = 2 so 
that 38 < n < 81 and the condition of case a is satisfied. 
In this case 81 - 38 = 5 x 8 + 3 so that a = 5 and /3 = 3 
and the construction Eq.lfTOli. i.e, [8 • 3] [> [14], gives rise to a 
stabilizer code [[38, 30, 3]], which is not optimal but the best 
code constructed so far The situation is similar for lengths 
n = 39, 82, 83. If n = 81 we have 

[81] - [2^] o [17], (12) 

an optimal code [[81, 73, 3]] obviously missing from the public 
code table. If n = 371 then m = 3 and 340 <n< 677 with 
the condition of case b satisfied. In this case 677 — 371 = 8 x 
38 + 2 so that a = 38 and (3 — 2 and by construction Eq. lfTTT i 
we have [371] = [8 • 26] > [2'^] [> [35] which is an optimal code 
[[371, 360, 3]] that saturates the quantum Hamming bound. 
Remarks For any given n > 38 we have either construction 
a or construction b. Generally: 

i) In the case of 8/„i + 3 < n < /,„+2 - 4 or /„i+2 + 

1 < n < 8/„i+i — 3 (unlabeled lengths in Table I) we have 
m = L^^^T^J ^'^'1 Theorem 2 gives rise to an optimal code 
because the quantum Hamming bound is saturated, i.e., its 
stabilizer has sh generators. 

ii) In the case of n = 8/™ or n = fm+2 with m > 2 
(lengths labeled by p in Table I) the construction of Theorem 

2 gives rise to a suboptimal code that is one logical qubit less 
than the optimal code given in Eq.© or Eq.(|7|i. 

iii) In the case of n = 8fm — {1, 2} or n = /rn+2 — {2, 3} 
with m > 2 (lengths labeled by u in Table I) Theorem 2 
gives rise to a so-far optimal code because in these cases its 
stabilizer has sh + 1 generators and there is no better code 
known so far However there may exist better codes with one 
more logical qubit. 

iv) In the case of ri = 8f„i + {1, 2} or n = fm+2 ^ 1 with 
m > 2 (lengths labeled by I in Table I) Theorem 2 gives rise 
also to a code with s^ + l stabilizers, i.e., one logical qubit less 



than that specified by the quantum Hamming bound. However 
it can be proved via linear programming bound that the codes 
constructed via Theorem 2 are optimal in these cases. □ 

IV. The linear programming bound 

In this section we shall work out analytically the LP bound 
for three families of lengths. For a stabilizer code [[n, fc, d]] 
we denote P as its projector and K = 2^ and s ^ n — k. The 
weight distributions Aj [14], [11] are defined by 

^» = ]^ E (» = 0,l,...,n). (13) 

where the summation is over all errors supported on i qubits. 
It is obvious that > 0, = 1, and J2i At = 2^* so that 
{^Ail2^Yi=\ be regarded as a probability distribution. For 
an arbitrary function f{x) we denote its average by 

1 

(14) 

■i=0 

In the following we shall formulate a subset of the linear pro- 
gramming bound for 1 -error correcting code, which serve our 
purpose perfectly. For a complete set of linear programming 
bound see Ref.[4], [12]. 

Linear Prograimning bound (Restricted set) If there exists 
a stabilizer code [[n, fc, 3]] then the following conditions hold 
true 

Ai = (3n-4a;), (15) 
^2 = i((4a; - 3?! + 1)2 -3n- 1), (16) 

5Z^2,>2^-\ (17) 

4=0 

Conditions Eq.dTSll and Eq.(fT6ll come from the error- 
correction conditions and condition Eq.lfTTIi comes from the 
fact that the even-supported subset of a stabilizer is a half or 
the whole stabilizer 

Theorem 3 If there exists a stabilizer code [[n, fc,3]], degen- 
erate or non-degenerate, for a length n equals fm+2 — 1 or 
8fm. + I or 8fm + 2 with fm = ^ and m = 1,2,..., 
then n — fc > Sff + 1, while the quantum Hamming bound 
for the 1 -error-correcting stabilizer codes of length n reads 
n ~ k > SH = [log2(37i + 1)]. 

Proof: Suppose that n = fm+2 ^ 1 with m > 1 it is 
obvious that ^ — 4™+^ — 1 is an integer and therefore the 
following function 

Ma;) = 16 (^x - ^) (^x - 1 - ^ j (18) 

is nonnegative for all integers. Specifically h{x) > for x = 
0,1,2, ... ,n. Furthermore we have 

/i(0) = 3?7.(3n + 4), (19) 
h{l) = 3??.(3n-4) > 2(3n + 4) (n > 3), (20) 
h{2) = (3?7.-4)(3n- 8) > 2(3n + 4) (n > 3). (21) 
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In this case the quantum Hamming bound reads sh = 2m+ 
4 (since 3n+4 — 2^™+'^) and we shall prove n — fc > s/j + 1 if 
there exists a stabilizer code [[n, fc,3]], which means Eqs.dTSl 
[TtT i must hold. As a result we have 

(/i(a;)) = 3n + 2Ai + 2^2. (22) 

Thus it follows from 

n 

2''{h{x)) = ^ > h{Q) + + h{2)A2 

> (3n + 4)(3n + 2Ai + 2A2) 

= (3n + 4)(/i(x)) (23) 

that 2* > 3n + 4, i.e., n — k > sh- We shall now prove that 
the equality can never happen. If the equality were true, i.e., 
2" = 3n + 4, then all the inequalities in Eq.(l23]l would become 
equalities which means that Ai ^ except i = 0,1,1 + 1 where 
I = 3n/4: since 1,1 + 1 are the only zeros of h{x). From 
conditions 1 + Ai + = 3n + 4 and 3n — 4j4;+i = 0, 

which comes from Eq.(fT5]l, we can solve Ai+i = in/A and 
Ai = Qn/A + 3. Noticing that / = 4'"+^ - 1 is odd it follows 
from inequality Eq.(fT7]i that 1 + 3n/4 > [in + 4)/2 which is 
impossible. Thus 2" > 3n + 4, i.e, n — k > sh + 1- 

The cases n ~ 8fm + {li 2} have been proved in [17] here 
we shall reproduce them for completeness. At first we suppose 
n = 8fm + 1 with m > 1. In this case sh = 2rn + 3 and we 
introduce a nonnegative function as 

= 16 (.T-^)'. (24) 

It is easy to check that as long as 77, > 5 

/(O) = (3n + l)2 > (3n + 5)(3?i-7) + 16, (25) 
/(I) = (371-3)2 > 4(3n + 5), (26) 
/(2) = (371-7)2 > 2(3n + 5) + 16. (27) 

If there exists a stabilizer code [[71, fc, 3]] then Eqs. dlSHlTt 
must hold. As a result of Eqs.( ll5m"6] ) we have 

(fix)) =3n+l + 4Ai+2A2. (28) 

As a result of Eq.dTTIi we have 

I6A0 + I6A2 + /(2«)^2. > 16 ^ > 8.2^ (29) 

i=2 i=0 

where we have used /(27) > 16 since '^"^^ , the unique zero 
of f{x), is an odd integer. Putting all these pieces together 

n 
i=0 

4/™ 

> /(0)+/(l)Ai+/(2)yl2+E/(2i)A2, 

> /(O) - 16 + /(l)Ai + (/(2)'-'l6)A2 + 8.2^ 

> (3?i + 5)(3n- 7 + 4^1 + 2/12) + 8.2'' 

= (3?i + 5)(/(x) -8) +8.2^ (30) 

in which the strict inequality comes from the /(O) term. As 
a result we have 2" > 3?i + 5 = 2"", taking into account of 
{f{x)) > 8. That is equivalent to saying n — fc > s/f + 1. 



Now we suppose ti = 8/,„ + 2 with 771 > 1 and in this case 
Sh = 2m + 3. Since ^2^2. is an integer the function defined 
as follows 

,(., = 16(.-2!i±l)(.-5i^) ,31, 

is nonnegative for integer x. It is obvious that g{0) > (3n + 
2)(3n-4) and g{i) > 2(37i + 2) fori = 1,2 as long as 71 > 5. 
If there exists a stabilizer code [[n, fc, 3]] then Eqs. (ll5Hl"7] i must 
hold, which leads to 

(<7(a;)} = 371-4 + 2^1 +2A2. (32) 

Thus 

r{9{x)) > 5(0) +5(lMi +.9(2)^2 

> (37i + 2)(37i-4 + 2^i + 2A2) 

= (37i + 2)(5(x)). (33) 

The strict inequality sign is due to the g{0) term. Since 
(5(2;)) > we have 2''' > 37i + 2 = 2'*", i.e., 71-fc > s^ + l. 

■ 

V. Special constructions 

In this section we shall prove Lemma 1 by constructing 
explicitly all the optimal non-degenerate codes with lengths 
71 < 37 except n = 6. Our main tool is a generalization of 
the pasting of stabilizer codes to a pasting of 2-error detecting 
blocks (2ed-block) defined as below. 

Definition 4 A 2-error detecting block [71, s]e is generated by 
a set of s multilocal Pauli operators acting on n qubits with e 
pairs being non-commuting that detects up to 2-qubit errors. 

Each non-degenerate stabilizer code [71, s] detect all 2-error 
and so they define 2ed-blocks [n, .s]o with all the generators 
being commuting. By shortening a pure code we generally ob- 
tain 2ed-blocks with some noncommuting pairs of generators. 
Some examples of 2ed-blocks are presented in Table V. 

2ed-blocks pasting: Given two 2ed-blocks [112,32] e2 
[iii,si]ej that are generated by (5*1 = X(7i2),S'2 = 
Z{n2), . . . ,Ss2) and {Ti,T2 . . . , T^^) respectively, then s = 
maxjsi, S2+2} generators as given in Table III is a 2-ed block 
[7ii+7i2,s]e with |ei— 62! < e < ei+e2- For convenience we 
shall denote by [iii, S2]ei 1> [^2, sijea the resulting 2ed-block. 

TABLE V 

Some examples of 2-error-detecting blocks. 





X 


/ XXX 


X X X X 






Z 


I z z z 


z z z z 






I 


X XYZ 


X Y Z I 






I 


Z Y Z X 


Y Z X I 






I2,4j2 13,4]2 


[4,4]i 




XXX 




X X X X X 


X X X X X X X 


z z z 




z z z z z 


z z z z z 


Z Z 


Z I z 




Y X Y X I 


Z I Z Y X 


Y X 


Z X Y 




I Z X Y I 


Z X Y I Z 


X Y 


Y Z X 




Z X I Y I 


Y Z X Z X 


I Y 






15,5|2 


[7,5]i 
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The 2ed-block given in Table III detects up to 2-qubits 
errors because firstly all the errors happening on the rii -block 
or ??2-block can be detected because [ni,si]ei and [712,52] 
are two pure codes of distance 3 and secondly two qubits 
errors happening on different blocks can be detected by the 
first two generators X{n2) (E) I{ni) and Z{n2) ^ I{ni). If 
two noncomuting generators are arranged in the same row 
the resulting generators will become commuting. As a result 
e can be zero when ei = 62 and all noncommuting pairs 
are carefully matched. In this case we obtain a pure 1 -error- 
correcting stabilizer code, since all 2-qubit errors can be 
detected. 

From the above arguments we see that although the 1 -qubit 
block, denoted as = {X,Z), detects only single qubit 
errors, it can be regarded as a 2ed-block because there is no 
2-qubit errors on a single qubit block. For example we have 
[2,4)2 = [l]i l> [l]i- As another example the perfect code 
[[5, 1, 3]] in Eq.Q can be regarded as the pasting of two 2ed- 
blocks [4,4]i [> [l]i. 

A 2ed-block fails to define a code because there are 
some pairs of noncommuting generators. By pasting two or 
more 2ed-blocks these noncommuting generators may become 
commuting and we thus obtain a 1 -error correcting stabilizer 
code. Our construction is therefore a kind of puncturing 
plus pasting. By puncturing some old stabilizer codes we 
obtain some 2ed-blocks that generally contains some pairs of 
noncommuting generators. By pasting with some other 2ed- 
blocks and carefully matching their noncommuting pairs we 
are able to produce some new stabilizer codes. To complete 
the constructions given in Table IV we have only to construct 
explicitly all the relevant 2ed-blocks. 

We consider the optimal code [2^] as in Table VI whose 
stabilizer is defined by the check matrix [RH^\Ac,RHc,] with 



/ 1 1 \ 
110 10 
10 
110 1 

V 1 1 ) 



/ 1 \ 

1 

10 

1 

\ 1 1 1 J 



(34) 



Obviously A5 is revertible and fixed-point free and R is in- 
vertible. By removing four coordinates [cs, cio, cig, 023] from 
this [2''] we obtain the 2ed-block [28,7]2 and by removing 
the first four coordinates [cq, ci, C2, C3] we obtain A 2ed-block 
[28, 7]i. By 2ed-blocks pasting with 2ed-blocks in Table V we 
obtain the pure optimal codes of lengths 30, 31, 33 and 35 in 
addition to a previously unknown optimal code 



[28,7]2C>[7,5]i>[l]i 



(35) 



whose stabilizer is explicitly given in Table VI. 

From three partitions of [2^] as shown in Table VII we 
can obtain a pure optimal code [10] as well as the unique 
optimal code [[6,0,4]] of distance 4 and four different 2ed- 
blocks. By pasting with the perfect 5-qubit code we obtain 
[15] = [10] O [5]. Also we obtain all the optimal pure codes 
of lengths from 11 to 14 as well as an optimal pure [7] = 
[6, 6]i [> [l]i. Finally the remaining 2ed-blocks appeared in 
Table IV are given in Table VIII. 



TABLE VII 

Three partitions of the optimal code [2'*]. 



[24] = [[16, 10, 3]] 



xxxxxxxxxx 
zzzzzzzzzz 

I X Y Z I I I X Y Z 

I Y Z X I I I Y Z X 

I I I I X Y Z X Z Y 

I I I I Y Z X Y X Z 



X X X X X X 

z z z z z z 

Y X Z Z Y X 
Z Y X X Z Y 
X Y Z X Y Z 

Y Z X Y Z X 



[10] = [[10,4,3]] 



[[6,0,4]] 



XXXXXXXXXX 

zzzzzzzzzz 

IXYZYXZZYX 
I YZXZYXXZY 
I I I I X Y Z X Y Z 
I I I I Y Z X Y Z X 



X X X X X X 

z z z z z z 

I I I X Y Z 

I I I Y Z X 

X Y Z X Z Y 

Y Z X Y X Z 



[10, 6] 1 



[6,6] 



XXXXXXXXXX xxxxxx 



zzzzzzzzzz 

I X Y Z I I I X Y Z 

I Y Z X I I I Y Z X 

I I I I X Y Z X Y Z 

I I I I Y Z X Y Z X 



z z z z z z 

Y Z X Z X Y 
Z X Y X Y Z 
X Y Z X Y Z 

Y Z X Y Z X 



[10,6]2 



[6,6]2 



TABLE VIII 
Further constructions of 2ed-blocks. 



[5,5]2 


[5, 5] 2 [5, 5] 2 


[3,5]2 


[7,5]i 


[5, 5] 2 [3, 5] 2 


[3, 5] 2 


7(5) 


X(5) y(5) 


Z{3) 


7(7) 


X(5) y(3) 


Z(3) 


7(5) 


y(5) z(5) 


X{3) 


7(7) 


Y(5) Z{3) 


X(3) 




[18, 7] 1 






[18, 7] 2 




[7,5]2 


[5,5]2 [5,5]2 


[3,5]2 


[7, 5] 2 


[7, 5] 2 [7, 5] 2 


[5,5]2 


7(7) 


X(5) Y{5) 


Z(3) 


7(7) 


X{7) Y{7) 


Z(5) 


7(7) 


Y{5) Z(5) 


X(3) 


7(7) 


Y{7) Z{7) 


X(5) 




[20, 7] 2 






[26, 7] 2 





VI. Discussions 

We prescribe a general construction of all the optimal 
stabilizer codes of distance 3 for lengths > 38 by past- 
ing known codes and a special construction of the optimal 
pur stabilizer codes of length 5 < n < 37 case by case 
by employing a generalization of the stabilizer pasting to 
noncommuting set of stabilizers, i.e., 2ed-blocks pasting. For 
three families of lengths we have worked out analytically the 
linear programming bound, which is strictly stronger than the 
quantum Hamming bound and ensures the optimality of our 
codes for these lengths. Except n = 6 all the optimal codes 
are pure. 

Apparently the construction given by Theorem 2 is not 
unique. Firstly there are different constructions for the optimal 
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TABLE VI 

The stabilizer for the optimal code [[36, 29, 3]] . 



[[32,25,3]] 



5 


10 


19 


28 





1 


2 


3 


4 


6 


1 


8 


9 


11 


12 


13 


14 


15 16 


17 


18 


20 


21 


22 


23 


24 


25 


26 


27 


29 


30 


31 
















X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


I 


I 


I I 


/ 


/ 


I 


I 


Z 


z 


z 


Z 


Z 


z 


z 


Z 


Z 


z 


z 


Z 


z 


z 


Z 


z 


z 


z z 


Z 


z 


Z 


Z 


z 


Z 


Z 


Z 


Z 


Z 


Z 


Z 


z 


I 


I 


I I 


/ 


/ 


I 


I 


z 


I 


X 


Y 


I 


z 


I 


z 


I 


I 


z 


I 


z 


z 


I 


z 


I 


Z Y 


X 


Y 


Y 


X 


Y 


X 


Y 


X 


Y 


X 


X 


Y 


X 


X 


X 


X X 


X 


X 


X 


I 


Y 


z 


Y 


Z 


I 


Y 


z 


X 


I 


z 


X 


I 


Y 


X 


I 


Y 


z 


X Z 


X 


I 


Z 


X 


I 


Y 


Z 


X 


I 


Y 


X 


I 


Y 


z 


z 


z z 


z 


z 


z 


I 


Z 


X 


Z 


X 


I 


Z 


I 


z 


I 


I 


z 


X 


Y 


Y 


X 


Y 


X 


Y I 


z 


I 


I 


Z 


I 


Z 


X 


Y 


X 


Y 


Y 


X 


Y 


Y 


z 


X z 


X 


I 


Y 


X 


I 


X 


X 


I 


I 


Z 


X 


Y 


Z 


Y 


X 


I 


Z 


Y 


z 


I 


Y 


X Z 


I 


Y 


I 


Z 


X 


Y 


Z 


I 


Y 


X 


Z 


X 


Y 


Z 


X 


Y I 


z 


X 


Y 


z 


Y 


Y 


Y 


Y 


I 


Z 


I 


Z 


X 


X 


Y 


Y 


X 


X 


z 


I 


Z 


I X 


Y 


X 


I 


Z 


I 


Z 


Z 


I 


Z 


I 


X 


Y 


X 


Z 


I 


Z Y 


X 


Y 


X 


I 




































[28,7]2 






























[7,5] 


1 









code [2™] [4]. Secondly there are other constructions such as 

[8 • (22"-i - ai)] > [8 • (22™-3 _ ^2)] ^ ^ _ 

...[>[8-(23-a„_i)]c>[17-/3] (36) 

or 

[8 • (2^™ - ai)] O [8 • (22'"-2 _ as)] C> . . . 

...o [8 -(24 -«,„_!)]> [37-/3] (37) 

where ai + 3 < 22(™-*+i)-i or 22('"-*+i) respectively and 
a = X^ta^ o^i- Poi" different choices of {ai} the resulting 
codes may be inequivalent. One problem arises as to the 
classification of all the optimal codes. Finally our approach 
can be easily adapted to nonbinary codes. 
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